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ABSTRACT 

High-resolution imaging has revealed an unusually high binary fraction 
amongst objects spanning the transition between the L dwarf and T dwarf spec- 
tral classes. In an attempt to reproduce and unravel the origins of this apparent 
binary excess, I present a series of Monte Carlo mass function and multiplicity 
simulations of field brown dwarfs in the vicinity of the Sun. These simulations 
are based on the solar metallicity brown dwarf evolutionary models and incorpo- 
rate empirical luminosity and absolute magnitude scales, measured multiplicity 
statistics and observed spectral templates in the construction and classification 
of composite binary spectra. In addition to providing predictions on the num- 
ber and surface density distributions of L and T dwarfs for volume-limited and 
magnitude- limited samples, these simulations successfully reproduce the observed 
binary fraction distribution assuming an intrinsic (resolved) binary fraction of 
11^3% (95% confidence interval), consistent with prior determinations. However, 
the true binary fraction may be as high as 40% if, as suggested by Liu et al., 
a significant fraction of L/T transition objects (~66%) are tightly-bound, unre- 
solved multiples. The simulations presented here demonstrate that the binary 
excess amongst L/T transition objects arises primarily from the flattening of 
the luminosity scale over these spectral types and is not inherently the result 
of selection effects incurred in current magnitude-limited imaging samples. In- 
deed, the existence of a binary excess can be seen as further evidence that brown 
dwarfs traverse the L/T transition rapidly, possibly driven by a nonequilibrium 
submergence of photospheric condensates. 

Subject headings: Galaxy: stellar content — methods: numerical — binaries: 
visual — stars: low mass, brown dwarfs stars: luminosity function, mass function 
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Introduction 



L dwarfs and T dwar fs are the two loAv est-luminosity spectral classes of low mass stars 
and brown dwarfs known fl Kirkpatrickll2005l . and references therein), span ning effective tem- 
peratures 700 < Tpff 5 , 2300 K and luminosities -5.7 < log Lf,oi/ Lq < -3.5 (iGolimowski et al. 



20041 : IVrba et al .1120041 ). Their photometric and spectral properties are remarkably distinct. 
L dwarfs typically exhibit red near-infrared colors {J — K 1.5 — 2.0), largely due to the 
presence of optically thick condensate dust in their photospheres. L dwarf spectra are char- 
acterized by weak bands of metal oxides (with the exception of CO), strong bands of metal 
hydrides (including FeH, CrH and H2O) and strong neutral metal lines, in particular al- 
kali species. Cooler T dwarfs, on the other hand, appear to have little or no condensate 
dust in their photospheres, exhibit prominent CH4 and H2O bands and pressure-broadened 
H2 absorption in their near-infrared spectra, and generally have blue near-infrared colors 
(J — K ^ — 0.5). There are ~450 L dwarfs and ~100 T dwarfs currently kno wnjll most iden- 



tified in wide-field imaging surveys such as the Two Micron All Sky Survey (ISkrutskie et al. 



2006 



1997 



hereafter 2MASS), the Deep Near-Infrared Survey o f the Southern Sk y (lEpchtein et al. 



hereafter DENIS) and the Sloan Digital Sky Survey (lYork et al.ll2000l . hereafter SDSS). 



The disparate properties of L dwarfs and T dwarfs have fostered attention on the tran- 
sition between these two spectral classes, a transition that has proven to be rather un- 
usual. It was recognized early on that the range of effective temperatures spanning the L/T 



transitio: 



is sm all, perhaps 100-300 K flKirkpatrick et all I2OOOI : IGolimowski et al.l 12004 



Vrba et al.l l2004j ). This is surprising, given the dramatic changes in spectral morphology 
that occur over this range. Parallax measurements of field sources further revealed an 
unexpected brightening of ~1 mag between the latest- t ype L dwarfs and earliest-type T 
dwarfs in the J-band (A^ ~ 1.2 /xm; iDahn et al. J2002h : iTinney. Burgasser. Sz Kirkpatrickl 



(I2OO3I ): IVrba et al.l (120041 )). This is also surprising, given the slow decline in bolometric lu- 
minosities over this range. Other unexpected trends across the L/T transition include the 
observed resurgence in gaseous FeH bands, a species that is expect e d to be depleted by the 
formation of condensates in the L dwarfs (IBurgasser et al.l l2002al : iMcLean et al.l |2003| ) ; a 
restrengthening of J-band K I lines after weakening considerably in the lates t-type L dwarfs 
(IBurgasser et al.l l2002l : iMcLean et al.l l2003l : ICushing. Rayner fc Vaccall2005l): and enhanced 
CO absorption at 4.7 /xni , possibly due to vertical upwelling (INoU. Geballe. &: Marleyl Il997l : 
Oppenheimer et al.l Il998l : iLeggett et al.l 120021 : IGolimowski et al.l |200J) . The current genera- 
tion of brown dwarf atmosphere models incorporating condensate clouds have been unable 



^ See |http : / / www . dwarf archives . org 

^The L/T transition is defined here as spanning spectral types L8-T4; cf. Leggett et al. ( 200Cll ). 
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to reproduce these trends (lAckerman &: Marlevll200ll : ICooper et al.ll2003l : iTsuji fc Nakajima 



2003 : Burrows. Sudarsky &: Hubenv 2006[). so dynamical (nonequilibrium) processes such 



as clou d fragmentation (jBurgasser et aLll2002al ). a change in condensate sedime ntation effi- 
ciency ( iKnapp et al.ll2004l ) and a global collapse of the condensate cloud layer (lTsujill2005l ) 
have been evoked as possible drivers of the L/T transition. However, persi stent empirical un- 
certa inties, such as the unknown ages an d surface gravities of field dwarfs (ITsuii fc Nakajima 
20031 ) and the role of unresolved binaries (IBurrows. Sudarsky fc Hubenyll2006l ) have obviated 
conclusive results. 

Recent high resolution imaging studies have helped to clarify our empirical understand- 
ing of the L/T transition. These studies have identified binaries whose components straddle 



the transition ( Cruz et al. 



2004 



McCaughrean et al 



2004 



Burgasser. Kirkpatrick fc Lowrance 



20051 : iBurgasser et al. 



2005 



2006a 



Liu et al 



12006 



Reid et al.l l2006bl ). two of which 



SDSS J102109.69-030420.1AB, a Tl + T5 pair (hereafter SDSS 1021-0304); and SDSS 
J153417.05+161546.1AB, a T1.5 + T5.5 pair (hereafter SDSS 1534+1615) — have been 
particularly revealing. In both systems, the later-type secondary is the brighter source at J- 
band, indicating that the brig htening trend observed i n the parallax d ata is truly an intrinsic 
feature of the L/T transition (jBurgasser et al.ll2006al : iLiu et al.ll2006l ). However, the bright- 
ening between the binary compo i ients, o f ord er ^0.2 mag, is si gnificantly less than previously 
inferred. Both IBurgasser et al.l (l2006a] ) and iLiu et al.l (120061 ) have hypothesized that prior 
results have been biased by unresolved multiple systems. Indeed, IBurgasser et al.l ( l2006al ) 
have demonstrated that the binary fraction of late-type L and early-type T dwarfs is as high 
as 40%, twice that of earlier-type L and later-type T dwarfs. This is likely a lower limi t given 
the probable existence of more tightly bound, unresolved systems. Indeed, iLiu et al.l (120061 ) 
hypothesize that most apparently single early-type T dwarfs may have unresolved compan- 
ions. Both studies ascribe this binary excess to the contamination of magnitude-limited 
samples by systems consisting of earlier-type a nd later-type component s, whose composite 
spectra mimic an L/T transition object. While IBurgasser et al.l ( l2006al ) demonstrated that 
the binary excess can be roughly reproduced assuming a rapid evolution of brown dwarfs 
across the transition, several additional factors — including the uncertain cooling rate and 
absolute magnitude scale of brown dwarfs across the transition, the underlying mass ratio 
distribution and frequency of brown dwarf binaries, and possible systematic effects arising 
in magnitude-limited samples — were not considered. Hence, the origin of the binary excess 
and its implications on the properties of the L/T transition remain unsettled questions. 

In this article, I present Monte Carlo mass function and binary spectral synthesis simu- 
lations aimed at examining the origin of the L/T transition binary excess and its dependence 
on various underlying factors. The empirical multiplicity dataset is described in detail in § 2. 
In § 3, the construction and implementation of the simulations are described, including real- 
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izations of the underlying mass function, binary frequency, mass ratio distribution, empirical 
luminosity and absolute magnitude trends, and the classification of unresolved binaries. Re- 
sults are presented in § 4, focusing on number density, surface density and binary fraction 
distributions for local volume-limited and magnitude-limited samples, and their dependen- 
cies on various input parameters. Results are discussed in § 5, including a breakdown of 
the origin of the binary excess and an examination of the intrinsic brown dwarf binary frac- 
tion based on the possible presence of overluminous, unresolved sources in current parallax 
samples. Conclusions are summarized in § 6. 



2. The Observed Binary Fraction 



The empirical multiplicity sample used here is similar to that constructed in lBurgasser et al 



( l2006al ). and is based on results from high resolution imaging studie s using the Hubble Space 
Telescope (HST) and ground-b a sed adaptive optic s instrumentation (IMartfn. Brandner fc Basri 



199a 



Reid et al 



20031: Gizis et al 



200j. These studies were chosen as they represent a relatively uniform set 



2001bl . I2OO6I : IBouv et all l2003l : iBurgasser et al.l 12003 



2006a 



Close et al. 



of observations in terms of angular resolution (limits of O'.'05-O'.'l) and imaging sensitivity. 
The selection criteria for sources imaged in these programs vary, but nearly all are composed 
of objects identified in color-selected, magnitude-limited surveys based primarily on 2MASS, 
DENIS and SDSS data. It is therefore assumed that, to first order, these data collectively 
constitute a magnitude-limited sample. Of course, there is no single magnitude limit nor 
filter band to which that limit corresponds that spans the entire sample. As I will show 
in § 4.2.3, this detail appears to be relatively unimportant in the resulting binary fraction 
distribution. 

For all of the sources observed in these imaging studies, spectral types were deter- 



mined from published optical classifications for L0-L8 dwarfs (on the iKirkpatrick et al. 



fll999h scheme) and ne ar- infrared data for L9-T8 dwarfs (on the iGeballe et al.l (|2002| ) and 
Burgasser et al.l (l2006al ) schemes). This division in classification stems partly from the wave- 
length regimes over which the L and T dwarf classes have been defined, and specifically the 
current absence of a robust near-infrared classification scheme for L dwarfs (see §3.5.2). Addi- 
tionally, few of the L dwarfs in these samples have reported near-infrared spectral data, while 
nearly^ll of them have been observed at optical wlvdengthsg Conversely, while all of the 



^Two exceptio ns to this are 2MASS J00452143+1634446 and DENIS J225210.7-173013, both from the 
Reid et al.l J2OO6I) HST study, fo r which only near-i nfrared spectral types have been reported (L3.5 and TO, 
respectively: Iwilson et all |2003h : iReid et all |2006l) ). 
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T dwarfs in these samples have been classified o n the iBurgasser et al 



schem e, few have been observed in the optical (IBurgasser et al 



2000 



2003a 



2006a ) near-infrared 



Liebert et al. 



2000l ). In cases where the near-infrared spectrum of a so urce indicated a spectral type of L9 or 



later, such as the L7.5/T0 SDSS J042348.57-041403.5 flGeballe et al.ll2002l : ICruz et al.l 12003 



her eafter SPSS 0423- 0414), the near-infrared classification was retained. Several sources in 
the IBouv et al.l (120031 ) HST sample required significant revisions to their reported classifica- 
tions (>3 spectral types) or were rejected given the absence of any published spectral data. 
Care was also taken to identify redundant sources between the samples. Finally, although th e 
L2 dwarf Kelu 1 was unreso lved in HST ob s ervat ions bvlMartm. Brandner fc Basril (119991) . 
more recent observations by iLiu fc LeggettI (120051 ) and iGelino. Kulkarni fc Stephens! (120061 ) 
have revealed this source to be binary, and is included as such in the sample. 

Table [T] lists binary fractions for the full sample in group s of integer subtypes, in cluding 
90% and 95% confidence intervals based on the formalism of IBurgasser et al.l (120031 ). While 
the overall sample is fairly large (32 binaries in 162 systems), significant uncertainties arise 
for individual subtypes due to small number statistics. The data were therefore binned into 
four spectral type groups, L0-L3.5, L4-L7.5, L8-T3.5, and T4-T8, chosen so as to minimize 
counting uncertainties while still sampling major spectral subgroups. There is a clear in- 
crease in the fraction of binaries from early-type L to the L/T transition, peaking at ef"^ = 
0-S3~^q'\1 for the latter grouping. Binaries are roughly twice as frequent amongst L/T tran- 
sition objects as compared to the earhest-type L dwarfs and latest-type T dwarfs. However, 
even with the prescribed binning, uncertainties due to small number statistics do not rule 
out a constant binary fraction of 0.20-0.24 at the 90% confidence level. Nevertheless, the 
simulations described below will demonstrate that a binary excess amongst L/T transition 
objects is an expected feature of the local brown dwarf population. 



Simulations 



3.1. Motivation and Underlying Assumptions 



Both IBurgasser et al.l (l2006al ) and lLiu et al.l (120061 ) have proposed that the frequency of 
binaries amongst L/T transition objects is a natural consequence of the properties of brown 
dwarfs in general. In order to properly investigate this hypothesis, it is therefore necessary 
to reproduce the empirical sample as observed. To do this, several basic assumptions were 
required. First, as nearly all L and T dwarfs known reside in the immediate vicinity of 
the Sun, it was assumed that the distribution of these objects is spatially isotropic. It was 
also assumed that the entire population arises from a single underlying mass function (MF) 
and age or birth rate distribution, as described in § 3.2. Third, only single sources and 
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binaries were modeled. Higher-orde r systems (triples, etc.) make up only 3% of all very 
low mass multiples currently known (IBurgasser et al.ll2006bl ). and are therefore a negligible 
population for the purposes of this study. Fourth, the modehng of binaries assumes a fixed 
overall binary fraction (e;,) and mass ratio distribution {P{q), where q = M1/M2). While 
there is evidence that the binary fraction of stars decreases across several stellar classes 



{ Duciuennoy fc Mayon Il99ll : iFischer fc Marcyl Il992l : IShatsky fc Tokoviniru l2002l : iBouy et al. 



2OO4J ). no trend has been established within the substellar regime. Fifth, while simulated 



populations were initially constructed in terms of space densities, they are re-sampled to 
reproduce a magnitude-limited sample in order to investigate systematic effects and make 
direct comparisons to the empirical data (see § 3.6). Finally, it is assumed that all binaries 
generated by the simulation are unresolved. This constraint arises from the observation 
that near ly all field L and T dwa rf binaries identified to date have apparent separations less 
than 1" (IBurgasser et al.l l2006bl lR. below the angular resolution limits of 2MASS, DENIS 
and SDSS. Hence, a given magnitude-hmited sample selected from these surveys consistents 
entirely of singles and unresolved multiples. 



3.2. Generating a Brown Dwarf Population: Mass Function and Birthrate 



The first step of the simulation, generation of a populat ion of low mass stars and brown 
dwarfs, was done using the Monte Carlo code described in Burgasser (|200^)i Briefly, this 
code generates a large set {N = 10^) of sources with randomly assigned masses and ages 
chosen from the MF and age distributions. These sources are then assigned luminosities 
and Tpf fS accordi n g to a set of evolutionary models, in this case from iBurrows et al.l (1200 ll ) 
and iBaraffe et al.l (120031 ). The distribution of luminosities can be compared to an observed 
luminosity function to constrain the underlying MF and birth rate (e.g., Burgasser 2001, 
2004; Allen et al. 2005). Here, the derived luminosity functions provide a starting point for 
examining the distribution of spectral types for single and multiple sources. 



The input parameters for the simulation are listed in Table [21 Power-law forms of the 



See http : //paperclip . as . cirizona. edu/$\sim$nsiegler/VLMjDinaries| Only one L dwarf binary 



wider than 1" is currently known, the L1.5+L4.5 pair 2MASS J15200224-4422419AB ([Kendall et al 



2006[ : [Burgasser et al.[ [2007[ ) . A handful o f wide binaries with late-type M dwarf p r imaries and L or 



T dwarf secondaries are also known (,Billeres et al. 2005 : Burgasser fc McElwainI 2006 : Biller et al. 20061 : 



McElwain fc Burgasseil[2006[ : [Reid et alJ[2006[ ) 



^Comparable brown dwarf popu la tion simula t ions h ave also been prese n ted in 



Chabrier (2002); Pirzkal et al. 


and Pinfield et al. ( 


2006h. For 



Reid et al 



|2005[ ): [Ryan et all |2005h : [Carson et al] ([2006[): [Peacon fc Hamblv 



For an alternative Bayesian approach, see Allen et al. ( 2005 ) 



1999 ): 



20061) : 
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MF, ^(M) = dN/dM oc M"" with a = 0.0 0.5, 1.0 and 1.5 were examined; as well as 



the log-normal distribution of IChabrierl (120021 . see also Miller & Scalo 1979). A fixed mass 
range of 0.01 to 0.1 Mq was assumed. The lo wer mass l imit is below the typical range 
of masses that comprise field L and T dwarfs (iBurgasserl |2004| ). while the upper limit is 
con strained by the maxi mum mass incorpo r ated in the evolutionary models. Note that 
the iBurrows et al.l (120011 ) and iBaraffe et al.l (120031 ) models use non-grey atmospheres and 
include condensate opacity as boundary conditions but not in the emergent spectral energy 
distributions (so-called "COND" models; lAllard et al.l (120011 )). Spectral models based on this 
assumptio n are largely consis t ent with the properties of mid-type M and mid- and late-type 
T dwarfs (IMarley et al.l Il996l : iTsuji. Ohnaka. fc Aokil Il996al : iTsujil 120051 ) ; but are generally 
inconsistent with those of late-type M and early-typ e L dwarfs in which condensate o pacity 
is important ( Jones fc Tsujilll997l : lAllard et al.ll200ll ). However. IChabrier et al.l ( 20001 ) found 
<10% difference in the evolution of luminosity and Teff between models that incorporate 
photospheric condensate dust and those that do not. This is a relatively sma ll deviation 



given current empirical constraints on the brown dwarf field luminosity function (ICruz et al. 



20031 ) and binary fraction distribution (Table [T]). Hence, evolutionary tracks that incorporate 
condensate opacity were not included in the simulations. 

The age distribution for the simulations assumes a constant birth rate over the past 
0.1-10 Gyr. B04 have demonstrated that significant deviations in the resulting luminos- 
ity function occur only for extreme forms of the birth rate, e.g., exponentially declining 
star formation rates or short star-forming bursts. These forms are generally inconsistent 
with the observed continuity in the Galactic MF between low- and high-mass stars; for- 
mation rates, kinematics and spatial distributions of planetary nebulae, white dwarfs and 
H II r egions; nucleosynthesis yields; and the metallicity and activ i ty distributions of G and K 
stars (IMiller fc Scaldll979l : ISoderblom. Duncan, fc Johnsonlll99ll : lBoissier fc Prantzoslll999l ). 
Simulations with different minimum ages (0.1, 0.5 and 1 Gyr) showed no significant differ- 
ences in the output distributions (see §§ 4.1.4 and 4.2.5). The metallicities of the simulated 
sources were assumed to be fixed at solar values, again a constraint of the available evolu- 
tionary models. 



Number de nsities frorn the si mulations were normalized to the field low- mass star (0.1- 
1.0 Mq) MF of lReid et al.l Jl999h . ^(M) = 0.35(g^)-i i3 pc^s y^-i^ implying a number 
density of 0.0037 pc~^ over the range 0.09-0.1 Mq. Because all of the distributions are scaled 
by this factor, adjustment to refined measurements of the low-mass stellar space density can 
be readily made. 
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3.3. Assigning Spectral Types: The Luminosity Scale 

Once a sample of sources was generated and luminosities derived, spectral types were 



from 


Golimowski et all J200^ 


(2004 


); 


Liu & Lemett (2005 



dwarf binaries Epsilon Indi Bab, Kelu lAB and SDSS J1021-0304AB, respectively — are 
shown in Figure [H and are restricted to sources with Mboi uncertainties <0.2 mag. As 
in the binary sample describ ed in ^ 2, spectral types for the 32 sources in this plot are 



based on the optical scheme of iKirkpatrick et al.l (119991) f or the L dwarfs (with the exception 



of the Kelu lAB comp onents; see iLiu fc Leggettl (120051 )) and the near- infrared scheme of 



Burgasser et al.l (l2006al ) for the T dwarfs. Figure [T] also delineates a sixth-order polyno- 
mial fit to these data, with coefficients listed in Table [31 This fit, which has a scatter of 
0.22 mag, does not include me asurements for the overluminous T dwarfs SDSS J125453.90- 



012247.4 (ILeggett et a 



( iBurgasser et al. 



20001 . hereafter SDSS 1254-0122) and 2MASS J05591914-1404488 
2000d . hereaft er 2MASS 0559-1404) . Desp ite being unresolved in high- 



resolution imaging observations (IBurgasser et al.ll2003l . l2006af ). these sources are sufficiently 
overlu minous for their spectra l types that they are likely to be cl osely-separated binary sys- 
tems (IGohmowski et all I2004J : IVrba et al.l l2004l : iLiu et all l2006l . see 5 5.2). Furthermore, 



were these two sources included in the fit, the resulting Mbo;/spectral type relation would be 
non-monotonic, obviating the ability to assign spectral types for a given luminosity. As it is, 
the Mboi relation exhibits a remarkable fiattening between types L8 and T4, decreasing by 
less than 0.5 mag over this spectral type range. This behavior is consistent with the small 
change in luminosity and Teff previously ascertained across the L/T transition. 

An alternative method for assigning spectral types is to use the simulated effective 
temperatures. There have been several studies on the Te// /spectral type scale for field L 
and T dwarf^; however, all of these rely on an assumed radius (or distributi on of radii) or 



spectral modeling, and deviations between these methods have been noted (ILeggett et al. 



200ll : ISmith et al.ll2003l ). Because bolometric luminosity measurements do not require any 
further theoretical assumptions, use of an M^o; /spectral type relation was favored. 



'^See Basri et al. (2000): KirkDatrick et al. (2000): Bureasser et al. (20021): Dahn et al. (2002): 


Leeeett et al.l (20021): Schweitzer et al. J2OO2I): 


Golimowski et al. 


(2004): Nakaiima et all ( 


2004h: Vrba et al. 


(2004): and Burgasser. Burrows & Kirkoatrick 


(2006). 
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3.4. Incorporating Binary Systems 



Binary systems were incorporated by selecting Nun = ^bN of the simulated sources 
as primaries, and creating secondaries with the same ages and masses assigned according 
to the mass ratio distribution. Binary fractions spanning < < 0.7 were examined, as 
well as two forms of P{q): a fiat distribution and an exponential distribution, P{q) oc g'", 
where F = 4. The latter form derives from statistical analyses of the observed distribution 
of very low m ass stellar and brown dwarf binaries, wh i ch exhibits a strong peak for equal- 
mass systems (IBurgasser et al.ll2006al : iReid et al.ll2006l : lAUenl 120071 ). Luminosities, effective 
temperatures and spectral types were assigned to the secondaries in the same manner as 
the primaries. To facilitate computation, secondaries were forced to have masses greater 
than the minimum mass sampled by the evolutionary models (0.01 M©). Again, for brown 
dwarfs with L and T dwarf luminosities, this constraint had a negligible effect on the final 
distribution of secondaries. 



3.5. Composite Spectral Types 



3.5.1. Generating Composite Spectra 

The binaries generated in these simulations have combined light spectra determined by 
the relative contributions of their primary and secondary spectral types. As t he absolute 



magnitude /spectral type relations of L and T dwarfs are dist inctly non-linear (jPahn et al. 



2OO2I : iTinney. Burgasser. &: Kirkpatridal2003l : IVrba et al.ll2004j ). and as the evolution of spec- 
tral features across this range is complex, the composite spectra were simulated explicity. 
This was done by combining near-infrared spectral templates selecte d from a large libr ary 
(88 sources) of L and T dwarfs observed with the SpeX spectrograph (IRayner et al.ll2003l ) in 
low-dispersion p rism mode. Details on the acquisition and reduction of the SpeX prism data 
are described in IBurgasser. Burrows fc Kirkpatrickl (|2006| . see also Gushing et al. 2003 and 
Vacca, Gushing & Rayner 2003). The full library will be presented in a future publication. 

The spectral templates chosen are listed in Table H] and their spectra shown in Fig- 
ure [21 Selection of T dwarf templates was straightfo rward, as they were sirn ply chosen from 
the primary and secondary spectral standards from IBurgasser et al.l (|2006al ). L dwarf near- 
infrared spectral standards, on the other hand, have not yet been defined. Furthermore, 
substantial variance amongst the near-infrared spectra of L dwarfs w i th identical optical 
classifications has been noted (iGeballe et al.l l2002l: iMcLean et al.l l2003l : iKnapp et al.l |2004J : 



Gushing. Rayner fc Vaccal 12005 



Burgasser et al.ll2007l ). These variations likely arise from 



differences in the abundance of condensates in the photospheres of L dwarfs, which is depen- 
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dent on and/or augmen t ed by surface gravity, metallicity and rotational effe cts ([Stephens 
2OO3I : iKnapp et allbooi burrows. Sudarskv fc Hubenvlbood : Ichiu et al-lbooei ). As it is be- 
yond tlie scope of this study to investigate these effects in detail, L dwarf spectral templates 
were con servatively chosen as sources with J — K ^, colors close to the mean of their spectral 
subtype (IKirkpatrick et al.ll2000l : IVrba et al.ll2004j ). and with no obvious spectral peculiarities 
or significant deviations between published optical and near- infrared classificat i ons. Three 
of the L dwarf templates are optical spectral standards in the IKirkpatrick et al.l (119991 ) clas- 
sification scheme. Three other L dwarf template s are binaries, but with components that 



have similar magnitudes and photo metric colors (iGolimowski et al.l |2004J : iReid et al.l 12006 



Stumpf. Brandner fc Henning 



20061 ). Finally, an L9 standard, 2MASS J03105986 +1648155 



( Kirkpatrick et al. I l2000h was included as a tie between the L and T dwarf classes (cf . iGeballe et al 



( 120021 )). As shown in Figure [21 these spectra vary smoothly in the evolution of the most 
dominant spectral features, including H2O, CO and CH4 absorption and overall spectral 
color. 

The spectral templates were then calibrated to absolute fluxes according to an adopted 
absolute magnitude scale. Three MKO0 M/^/spectral type relations and two Mj/spectral 
type relations were examined. The primary relation used in the simulat ions was derived by 



fitting an eighth-order polynomial to Mk datslj based o n photometry (jGeballe et al 



Leggett et al 



2003 



2002 



Vrba et al. 



2002 



Knapp et al.ll2004f) and parallax data (jPahn et al. iTinney. Burgasser. fc Kirkpatr 



20041 ) for 26 unresolved field L and T dwarfs with combined uncertain- 

ute magnitudes for the binaries Epsilon Indi Bab 



ties <0.2 mag; and comp onent abso 



(iMcCaughrean et al.ll2004r). Kelu lAB (iLiu fc Leggettll2005h . SDSS J0423-0414AB and SDSS 



J1021-0304AB (IBurgasser et al.ll2006al ). Figure [3] displays these data and the polynomial fit, 
which has a dispersion of 0.26 mag (coefficients are listed in Table [3]). Again, the overlumi- 
nous sources 2MASS 0559-1404 and SDSS 1254-0122 were excluded from the fit. The other 



absolute magnitude/spectral type relations were taken from iLiu et al.l (120061 ) . based on fits 
to the absolute magnitudes of field sources excluding either known ("L06") or known and 
possible binaries ("L06p"; see § 5.2). Figure [3] compares the three Mft-/spectral type rela- 
tions. The most sig nificant di f ferenc es between these lie in the L8 to T5 range, over which 
the two relations of iLiu et al.l (120061 ) differ by ~1 mag. The polynomial fit derived here is 
intermediate between these two relations. T he Mj/spectral type relations considered here 
are illustrated in Figure 4 of iLiu et al.l (120061 ). 



Once the template spectra were scaled to their respective absolute magnitudes following 



''Mauna Kea Observatory system; ISimons k, Tokunagal (|2002l ): iTokunaga. Simons fc Vaccal (|2002l ). 



See http : //www . j ach . hawaii . edu/$\sim$skl/LTdata. html 
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standard techniques (cf. ICushing. Rayner fc Vaccal (120051 )). composite spectra were created 
by adding the templates that bracketed the assigned spectral types of the binary components, 
weighted accordinglj]^. Examples of composite spectra for various combinations of primary 
and secondary spectral types are illustrated in Figure HI 



3.5.2. Classifying the Composite Spectra 



For binaries with large differences in their component spectral types, composite spec- 
tra can differ substantially in their spectral characteristics from the underlying standard 
sequence. Examples include the appearance of CH4 absorption at iJ-band a nd not K- 
band, as observed in the L 6: + T4: binary 2MASS J05185995-2828372AB JCruz et al. 



2004 : iBurgasser et al.l l2006al . cf. Figure Hj); or strong deviations in individual H2O absorp- 
tion bands. As tens to hundreds of thousands of composite spectra were generated in each 
simulation, an automated classification scheme was required After experimenting with var- 
ious techniques, it was determined that classification by spectral indices provided the most 
efficient and accurate method. 



For T dwarfs, classification by indices is straightforward, as IBurgasser et al.l (l2006al ) 
have defined five near-infrared indices (H2O-J, CH4-J, H2O-H, CH4-H and CH4-K) which 
track monotonically with T spectral type. For the L dwarfs, several near-infrared indices 



have been defined and compared agains t optical types (IReid et al.l l2001al : iTesti et al.l 12001 



Geballe et al.ll2002l : iMcLean et al.ll2003l ). the most successful of which are tied to the H2O 
absorption bands. However, a robust classification scheme for L dwarfs in the n ear-infrared 



is sti ll lacking, which can lead to problems parti cularly at the L/T trans ition (jChiu et al. 



20061 ) . I therefore examined the behavior of the IBurgasser et al.l (l2006al ) indices for both 
L and T dwarf spectra from the SpeX prism library, as shown in Figure [51 Following the 
convention of the empirical binary sample, spectral types for L0-L8 dwarfs are based on 
optical data while those for L9-T8 dwarfs are based on near-infrared data. There are strong 
correlations for three of these indices — the two H2O indices and CH4-K — that are roughly 
monotonic across full spectral type range. CH4-J and CH4-H are strongly correlated only in 
the T dwarf regime. Fourth-order polynomial fits to the indices were made where correlations 
were strongest, and coefficients are listed in Table [31 

Spectral types for the composite spectra were determined from these spectal index 
relations, as follows. First a spectral type was estimated from the H2O-J, H2O-H and CH4- 
K indices. If these indices indicated a T spectral type, or if CH4 absorption is present at 



For example, an L1.3+T4.7 composite spectrum was modeled as 0.7xLl + 0.3xL2 + 0.3xT4 + 0.7xT5. 
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iJ-band (indicated by CH4-H < 0.97), subtypes based on the CH4-J and CH4-H indices 
were also computed. The individual index spectral types were then averaged to derive an 
overall classification for the composite. Figure H] lists the classifications derived for the binary 
combinations shown there. 

To determine the robustness of this method, the L and T dwarf SpeX spectra used to 
derive the index relations were reclassified by the above technique and compared to their 
original classifications. Figure E] illustrates this comparison and typical deviations. Overall, 
spectral types derived from the indices agree with published values to within 0.6 subtypes, 
although the scatter is greater amongst L0-L8 dwarfs (0.9 subtypes) than L9-T8 dwarfs 
(0.3 subtypes) or L8-T4 dwarfs (0.4 subtypes). These deviations likely arise from the same 
photospheric condensate var iations that hav e inhibited the identification of near-infrared 



L dwarf spectral standards (IStephend l2003l ). Nevertheless, Figure E] indicates an overall 
accuracy of 0.5-1.0 subtypes can be attained for both L dwarfs and T dwarfs using a common 
set of spectral indices. 



3.6. Constructing Magnitude-limited Samples 

Once the parameters for all single and composite systems were derived, the simulated 
population was resampled into integer spectral type bins to derive distributions in number 
density and binary fraction. These distributions were normalized according to the assumed 
space density (§ 3.2), yielding a "volume-limited" sample. The resulting number densi- 
ties were then converted to surface densities for a magnitude-limited sample by computing 
effective volumes for each system, 

V.ff^i = 1ioO-6('"h™-A^>)+3 (1) 



( ISchmidtl Il968l ). where mum is the adopted limiting magnitude and Mi the absolute mag- 
nitude of the system (i.e., the absolute magnitude of a single source or the combined light 
magnitude for a binary system). For all simulations, m^j^ = 16 was adopted as a proxy. 
Surface densities to this magnitude limit were then computed as the sum of effective volumes 
over a given spectral type range, 

ieSpT 

S(SpT) = J2 2.424xlO-Ve//,, deg'^ (2) 

i 

where the numerical factor provides units of deg~^. 
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Results 



A total of 58 separate simulations were run to sample the various input parameters listed 
in Table [21 From each simulation, volume-limited space densities, magnitude-limited surface 
densities, and binary fraction distributions (for both volume-limited and magnitude-limited 
cases) as a function of primary, secondary and systemic (singles plus composite binaries) 
spectral type were constructed. I n the comparisons that follow, the baseline simulation is 
based on the iBaraffe et al.l (120031 ) evolutionary models, a power-law mass function with a 
= 0.5, eb = 0.1, the exponential mass ratio distribution, and the M/^/spectral type relation 
derived in 5 3.5.1. 



4.1. Number and Surface Density Distributions 

While the focus of this study is the binary fraction distribution of L and T dwarfs, several 
clues into its origin can be gleaned from examining the number density (for the volume- 
limited case) and surface density (for the magnitude-limited case) distributions and their 
dependencies on various underlying parameters. Results from representative simulations are 
provided in Tables O and El 



Distribution of Single Sources and Dependencies on the Underlying Mass Function 



Figure [7] plots the distribution of singles sources (i.e., for = 0) as a function of 
spectral type for volume-limited and magnitude-limited samples and the five mass functions 
examined. The overall shape of the volume-limited distributions is similar for all of the mass 
functions, with a plateaio amongst early- and mid-type L dwarfs, followed by a decline in 
number density from mid-type L to a minimum around T0-T2 (a factor of 4-5 decrease); 
then, a rapid rise toward the latest-type T dwarfs. The shape of thes e distribut i ons is 
simi lar in form t o the simulated field brown dwarf luminosity functions of iBurgasserl (120041 ) 
and I Allen et al.l (120051 ) . and arises from the combination of long-lived stellar-mass early-type 
L dwarfs (which accumulate over the 10 Gyr of the simulation), the rapid cooling rate of 
substellar-mass L dwarfs (reducing the number of late-type L dwarfs) and the slow cooling 
rate of late-type T dwarfs (increasing their numbers). The dip at T0-T2 is enhanced by the 
flattening in the M^o; /spectral type relation over this range, which spreads the luminosity 



^"^The slight rise in number densities from LO to L2 is likely an artifact of the maximum mass (0.1 M©) 
used in the simulations. 
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function thinly at these spectral types. In addition, mass functions with a > preferentially 
produce lower-mass brown dwarfs and hence contribute to the rise in number density amongst 
the lowest-luminosity brown dwarfs. Overall number densities increase for steeper power-law 
mass functions, as expected, ranging over 15-60% from types LO to T8 for Aa = 0.5. The 
log-normal mass function produces a distribution intermediate between the a = 0.5 and 1.0 
distributions. The number density of L0-L9 (T0-T8) dwarfs in these simulations range over 
0.006-0.011 pc-3 (0.009-0.03 pc'^) for a = to 1.5; i.e., there are 50-250% more T dwarfs 
than L dwarfs in a given volume for the mass functions examined. The number density of 
L/T transition objects, L8-T4, range over 0.0018-0.004 pc~^ for a = to 1.5, and are thus 
comparatively rare. 

In the magnitude-limited case, the number density decline from mid-type L to early- 
type T combined with the decrease in absolute i^-band brightness over this range results in 
a steep decline in surface densities, by a factor of ~ 120-200 from LO to Tl depending on the 
underlying MF. Surface densities then increase by ~60-70% over types T1-T4, driven by 
both the increase in number densities and the flattening of the Mi^/spectral type relation 
over this spectral type range. Beyond T4, surface densities again drop as the Mi^-Zspectral 
type relation steepens. The intrinsic brightness of L dwarfs implies that there are 50-80 
times more of these objects than T dwarfs in a given magnitude-limited sample, in stark 
contrast to the ratio of number densities]^ Again, steeper mass functions result in more L 
and T dwarfs, by roughly the same factor as the volume-limited case. 



4.1.2. Distributions of Primaries, Secondaries and Systems 

Figure [H] compares the spectral type distributions of primaries, secondaries and systems 
(singles plus binaries) in both volume-limted and magnitude-limited cases. The only clear 
difference seen between the primary and systemic number density distributions in the volume- 
limited case is a 10-15% increase in the number of T0-T2 dwarf systems. This deviation 
increases for higher binary fractions, and while relatively small nevertheless has a significant 
impact on the binary fraction distribution. For the magnitude-limited distribution, the 
contribution of binaries in the surface densities of T0-T2 dwarfs increases substantially, by 
~60% for eb = 0.1, implying that 40% of these objects are binary. This is consistent with 
the observed binary fraction for these spectral types (Tabled]). Increasing the underlying 



^^The current tally of known L dwarfs and T dwarfs, a ratio of ^4.5:1, largely reflects preferential efforts 
toward searches of later-type sources. If both classes were investigated to similar (JsT-band) magnitude limits, 
perhaps 5000-8000 L dwarfs would be known. 
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binary fraction raises the entire surface density distribution in the magnitude-hmited case, as 
the over luminous, unresolved binaries are sampled to larger distances and greater volumes. 
However, the largest increase is consistently seen in the surface densities of T0-T2 dwarfs, 
reaching a factor of 4 increase for e?, = 0.5. 

The shape of the number density distribution of secondaries in Figure [8] closely follows 
that of the primaries and composites, but is scaled roughly by the underlying binary fraction. 
There is a slight increase in the relative numbers of T dwarf secondaries, an artifact of the 
constraint that secondaries have equal or lower masses than their primaries. Since M dwarf 
primaries are not modeled here, there is an artificial deficit of L dwarf secondaries. This 
bias is greatly enhanced in the magnitude-limited case. Interestingly, the number of T dwarf 
secondaries is substantially larger than the number of primaries in this case, by a factor 
of 3 or more for spectral types T4 and later and e^, = 0.1. The vast majority of these 
secondaries are low mass companions to L dwarfs, as their relative infrequency (lying on the 
tail of the power-law mass ratio distribution) is outweighed by the much larger number of 
brighter L dwarf primaries in a magnitude-limited sample. A significant fraction of these 
systems consistent of old, low-mass stellar (L dwarf) -|- brown dwarf (T dwarf) pairs whose 
components straddle the hydrogen burning limit. These results indicate that, for a given 
magnitude limited sample, many more T dwarfs may be found as companions to higher mass 
brown dwarfs and/or stars than as isolated objectsjlj although uncovering these companions 
may be nontrivial. 



4- 1.3. Dependencies on the Underlying Absolute Magnitude Relation 



The spectral type distribution of a magnitude-limited survey is naturally influenced by 
the underlying absolute-magnitude relation. As illustrated in Figure [9l the most significant 
differences found in these simulations are largely constrained to L/T transition objects, 
where the three absolute magnitude/spectral type relations considered here vary appreciably. 
Most pertin ent are th e differ ences between simulations based on the two M^/ spectral type 
relations of iLiu et al.l (120061 ). which differ by up to a factor of 2-5 over the L8-T5 range. 
This is consistent with the ~1 mag difference between these relations, implying a factor 



^^There is an indication of this trend when one considers that of aU known T dwarfs with dis- 
tance measurements within 8 pc, two are is olated field objec ts (2MASS J04151954-0935066 and 2MASS 
J09373487+2931409: lBurgasser et al.l |2002h : IW3a et al.1 |2004l )l. w hile five are coinpanions to nearby stars 
fEpsilo n Indi Bab; Gl i ese 22 9 B, Gliese 570D and SCR 18 4 5-6357B:lNakaiima etld] |l995l) : iBurgasser eTal 
|2000b[ ): IScholz et al] |2003h : iMcCaughrean et all |2004l ): iBiller et all |2006[) ). However, selection effects, 
including concentrated efforts to find low mass companions to nearby stars, may be responsible. 
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of ~4 increase in volume sampled. The large differences in the surface densities imply 
that a complete magnitude-limited sample of L/T transition objects can provide a robust 
constraint on the underlying absolute magnitude relation. Similar results are also found 
between different Mj/spectral type relations. Figure [HI also illustrates the differences between 
J-band limited and i^-band limited samples. Surface density distributions in the former case 
show considerably more structure, with more pronounced local minima and maxima at types 
TO and T5, respectively. This reflects the color shift between L dwarfs and T dwarfs (from 
red to blue J — K colors) , the sharp e r turn over in Mj/spectral type relations across the L/T 
transition (cf. Figure 4 in lLiu et al.l (120061 ) ) and the increase in number densities beyond TO 
(Figure E]). 



4-1 .4- Weakly Correlated Parameters 

Several factors considered in the simulations contributed minimally to the overall shape 
of the number density and surface density distributions. These include the different minimum 
ages for the field population examined, spanning 0.01 to 1 Gyr, which only served t o decrease 



the nu mber of early-type L dwarfs by <15%. This is consistent with the results of iBurgasser 



( 120041 ). which demonstrated that only extreme modulations of the birthrate change the 
luminosity function of brown dwarfs appreciably. The two evolutionary models examined 
also produced similar spectral type distributions, with small deviations amongst early-type 
L dwarfs arising from different hydrogen-burning minimum masses derived by these models 
(0.075 Mq for Burrows et al. 1997 versus 0.072 M© for Baraffe et al. 2003). Constant and 
exponential mass ratio distributions produced nearly identical spectral type distributions for 
et, = 0.1, although 30% differences in the surface densities of T0-T4 systems are seen for 
et, = 0.5. Just as an increase in the underlying binary fraction increases the total effective 
volume sampled for a given subclass, the exponential mass ratio distribution, which favors 
equal-mass/equal-brightness binaries, samples a larger effective volume than a flat mass ratio 
distribution. However, this effect is considerably smaller than the increase in number and 
surface densities (particularly amongst L/T transition objects) incurred by an increase in 
the underlying binary fraction. While all of these differences between the simulations are 
statistically significant, they are considerably smaller than cur rent empirical uncertainties 



for L and T dwarf number and surface density measurements (jBurgasseii l200ll : ICruz et al. 



20031 ). and can be considered negligible for the purposes of this study. 
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4.2. Binary Fraction Distributions 

4-2.1. Reproducing the Binary Excess 

Binary fractions as a function of spectral type were computed by ratioing the number 
of binaries witliin a given systemic spectral type range to all simulated sources within that 
same range. Results from representative simulations for both volume-limited and magnitude- 
limited cases are given in Tables [7] and El respectively. Figure [10] displays the binary fraction 
distribution as a function of systemic spectral type for the baseline simulations. For both 
volume-limited and magnitude-limited samples, there is a considerable degree of structure in 
these distributions, deviating significantly from the underlying fiat form. The largest jump 
is seen amongst T0-T2 dwarfs, over which binary fractions increase by a factor of two in 
the volume-limited case. The remaining low- level structure, of order 30-50%, likely arises 
from composite spectral classification errorJ^ or, in the case of LO-Ll dwarfs, endpoint 
effects. However, these are significantly smaller than the strong peak in the binary fraction 
distribution at types T0-T2. 

In the magnitude-limited case, all binary fractions are consistently higher, a well- 
understood bias incurred by the larger volume sampled by overluminous, unresolved pairs0 
For the most part, this bias merely serves to enhance the structure seen in the volume-limited 
distribution and, most importantly, the fraction of binaries amongst T0-T2 dwarfs, which 
rises to 0.43. This result was anticipated by the increase in the surface density distribution 
between systems and singles as noted in § 4.1.2 (Figure [S]). 

That the fraction of binaries amongst T0-T2 dwarfs is higher in a volume-limited sample 
suggests that the intrinsic binary fraction of early-type T dwarfs is higher. However, it is 
important to make a distinction between the binary fraction distribution of systems and 
that of primaries. The latter exhibit a fiat distribution in both the volume-limited and 
magnitude-limited cases, consistent with the assumed fixed underlying binary fraction. A 
more appropriate conclusion is that early-type T dwarf systems are more frequently binary, 
when such systems are unresolved. 



^•^Simulations with alternate spectral templates were examined to ascertain whether this significantly 
influenced the binary fraction distribution results. To the limits of our spectral sample, only minimal 
changes in the distribution were found (of order 10%). Note that as dips are generally followed or preceded 
by peaks, it is likely that such systematic effects are averaged out in the coarser spectral binning used for 
comparisons to the empirical data. 

Following the formalism of iBurgasser et al. I (|2003l) . a mass ratio distribution of the form P{q) oc q"^ and 
intrinsic binary fraction of 0.1 implies an apparent binary fraction of 0.19, a factor of roughly two increase. 
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4-2.2. Comparison to Empirical Data and Constraints on the Underlying Binary Fraction 



Figure [TT] compares the magnitude-limited binary fraction distribution as a function of 
systemic spectral type to the binned empirical data of Table [H The agreement between 
these for the baseline parameters is quite remarkable. The steady rise in the binary fraction 
between early-type L dwarfs and L/T transition objects, the peak of the distribution in the 
latter case, and the subsequent drop in the binary fraction for the latest-type T dwarfs are 
all faithfully reproduced in both form and magnitude. 

Figure [11] also illustrates how the binary fraction distribution varies with the assumed 
underlying binary fraction, eb- It is seen that while increasing increases the scale of the 
distribution, it has little effect on its shape; in particular, the relative number of binaries 
between L/T transition objects and early-type L or later-type T dwarfs. 

Comparison between these distributions and the empirical data enable a constraint on e^. 
By linearly interpolating the distributions as a function of e^, and minimizing the uncertainty- 
weighted deviations between the data and simulation results, a best fit e^^* = 0.11 was 
found, with ranges of 0.09-0.15 (0.08-0.17) acceptable at the 90% (95%) confidence levels. 
These values are fully consistent with earlier estimates of t he volume-limited , resolved binary 



fract i on of L and T dwarfs, al so ranging over 0.09-0.15 (IBouy et al.l l2003l : iBurgasser et al. 



2003. 


2006a: 


Reid et al. 


2006) 



\t is important to remember that e^*^* is the underlying 
resolved binary fraction, having been constrained by observations of resolved systems. The 
difference between this and the intrinsic binary fraction of brown dwarfs is discussed further 
in 5 5.2. 



4.2.3. Dependencies on the Mass Ratio Distribution 

A somewhat larger modulation of the apparent binary fraction distribution is seen when 
different mass ratio distributions are considered, as illustrated in Figure [131 A flat mass 
ratio distribution results in a 20-30% decline in the binary fraction distribution across the 
full spectral type range, in both the binned and unbinned cases. Again, this decline is 
consistent with the reduction in effective volume sampled by the reduced fraction of equal- 
mass/equal-brightness systems for a magnitude-limited survey. Figure [12] illustrates that 
a simulation using a fiat mass ratio distribution and = 0.14 produces a binary fraction 
distribution that differs by less than 10% from the baseline simulation. Indeed, the empirical 
data constrain a best fit e^*^** = 0.14 and acceptable ranges of 0.12-0.19 (0.11-0.21) at the 90% 
(95%) confidence limits in the case of a fiat mass ratio distribution, somewhat higher than 
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that derived using an exponential distribution^ This degeneracy between the mass ratio 
distribution and the underlying binary fraction implies that the former cannot be robustly 
constrained by the binary fraction distribution without an independent determination of e^. 



4.2.4- Dependencies on the Absolute Magnitude/ Spectral Type Relation 

Figure [12] illustrates how the binary fraction distribution varies for different underlying 
absolute magnitude distributions. Despite the substantial effect that the absolute magnitude 
relation has on the surface density distribution of L8-T5 dwarfs (Figure [9]), the impact on the 
binary fraction distribution is considerably smaller. Between the L06 and L06p Mj^/spectral 
type relations, binary fraction distributions differ by less than 20% for individual subtypes, 
and by less than 10% when the spectral types are binned according to the empirical data. 
There is a comparably small d ifference between binary fraction distributions generated by the 



Mj/spectral type relations of iLiu et al.l (120061 ) . and also between J- and K-haxid magnitude- 



limited samples. These comparisons indicate that the binary fraction distribution, and more 
importantly the binary excess at the L/T transition, are relatively insensitive to the details 
of how a magnitude-limited sample is constructed. This justifies the incorporation of several 
high-resolution imaging samples into a single magnitude-limited probe of the binary fraction 
distribution. 

A more important conclusion to draw from these comparisons, however, is that the 
observed binary excess does not appear to be an artifact of magnitude-limited imaging 
samples, since the parameters that most influence the surface densities of L and T dwarfs 
(the absolute magnitude scale and limiting filter band) have little bearing on the binary 
fraction distribution. Indeed, as already pointed out in § 4.2.1, a magnitude-limited sample 
simply amplifies the excess present in the underlying population. 



4.2.5. Weakly Correlated Parameters 

As with the number and surface density distributions, several of the parameters explored 
in these simulations had minimal infiuence on the binary fraction distribution. Neither the 
form of the mass function, choice of evolutionary model or minimum age of the population 



^^FoUowing the formalism of iBurgasser et al.l (j2003[ ). a flat mass ratio distribution and an inherent binary 
fraction of 0.1 implies an apparent binary fraction of 0.14, a decrease of 26% as compared to the exponential 
mass ratio distribution. 
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had more than a 5% effect on the structure or magnitude of the distribution. Out of all of 
the parameters varied, only e?, significantly influences the binary fraction distribution and 
only by scaling the entire distribution up or down. The structure of the distribution — i.e., 
the excess of binaries amongst L/T transition objects — must arise from some other aspect 
of the simulation. As discussed below, it appears that the luminosity scale is exclusively 
responsible for this excess. 



5. Discussion 



5.1. On the Origin of the L/T Transition Binary Excess 

The simulations presented here reproduce remarkably the relatively hig h binary fraction 



of L/T trans ition brown dwar fs. Yet how does this binary excess arise? Both lBurgasser et al. 



(l2006al ) and iLiu et al.l (120061 ) proposed that the similarity of the composite spectra of L + 
T dwarf binaries to early-type T dwarf spectral morphologies, and the rapid evolution of 
single brown dwarfs through the L/T transition, may be the underlying causes. A detailed 
examination of the simulations presented here support these hypotheses. 

First, the composition of L/T transition binaries in the simulations is quite distinct as 
compared to other spectral subclasses, as illustrated in Figure [TH Most binaries have com- 
posite types similar to their primaries, reflecting both the preference for equal- mass/equal- 
luminosity systems in the exponential mass ratio distribution and the dominance of the 
primary in the combined light spectrum of systems with low-mass companions. Yet binaries 
classified as L9-T3 dwarfs consistently have earlier-type primaries, up to ~0.6 subclasses on 
average, while 15-20% of T0-T2 systems have primaries at least one subclass earlier. This is 
consistent with studies of known L/T tra nsition binaries, 70% of which have been shown or 
inferred to be L dwarf plus T dwarf pairs (IBurgasser et al.ll2006al ). Such pairs readily mimic 
the spectral appearance of an early-type T dwarf. The flattening of the absolute magnitude 
scale across the L /T transition allows the T dwarf secondaries of these systems to contribute 
significantly to the overall spectral flux at near infrared wavelengths. Thus, a wide variety 
of component combinations result in a composite spectral energy distribution similar to an 
early-type T dwarf (cf. Figure H]). 

Second, a rapid evolution of brown dwarfs across the L/T transition can be deduced 
fr om the flattening of t he luminosity scale over this range. Based on the evolutionary models 
of iBaraffe et al.l (120031 ) and the Mfto^/spectral type relation shown in Figure (H a 0.05 M© 
(0.03 Mq) brown dwarf dims from type L8 to type T3 in ~600 Myr (~100 Myr), a period 
spanning 30% (20%) of its lifetime at that stage. In comparison, a brown dwarf of similar 
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mass spends 2-3 times longer cooling from LO to L8, and a full 7 Gyr (3 Gyr) cooling 
thereafter to the end of the T spectral class. Since single brown dwarfs spend relatively 
little time at the L/T transition, they are much rarer than earlier- type and later- type brown 
dwarfs in a given volume (Figure [7]). Yet, combinations of these other spectral types, even 
for small binary fractions, can be comparable in number to single L/T transition objects, 
resulting in the perceived binary excess. 

These two explanations for the high rate of binaries amongst L/T transition objects are 
in fact related to a single underlying cause: the flattening of the luminosity scale. The small 
decline in Mboi across the L/T transition directly translates into small changes in the absolute 
magnitude scale over this range, and is the root cause of the rarity of single L/T transition 
objects. Consider then the reverse argument: The observed excess of binaries at the L/T 
transition is further evidence of a flattening in the luminosity scale. Indeed, the agreement 
between the shape of the simulated binary fraction distribution and empirical data, which is 
largely independent of all other parameters besides the Mbo^ /spectral type relation, suggests 
that the luminosity scale of brown dwarfs across the L/T transition is indeed quite flat, 
and that brown dwarfs must traverse this transition over a relatively short period. This 
lends support to suggestions that the complete removal of photospheric condensate dust 
between late-type L and mid-type T, which appears to be largely responsible for the changes 
in spectral morphology across this transition, is relatively rapid and may require dynamic, 
nonequilibrium atmospheric processes. 



5.2. Absolute Magnitudes of Field L/T Dwarfs and the Intrinsic Brown Dwarf 

Binary Fraction 



In § 4.2.2, a constraint on the underlying resolved binary fraction, ejj*^**, was made by com- 
paring simulated binary fraction distributions to high-resolution imaging results. However, 
imaging cannot identify binaries that have very small separations or systems observed when 



(1999 


); 


Liu & Legeett 


(2005 


); 



tight brown dwarf binaries in the field is likely, as the separation distributi on of resolved sys- 



tems peaks at o r near the resolution lim it o f current imaging s urvey s ( iBouy et al.l 12003 



Gizis et all I2OO3I : ISurgasser eraP l2006al jbll | Maxted fc Jeffries! (120051). analyzing the re- 
sults of the first searches for spectroscopic brown dwarf bin aries (IBasri fc Martini Il999a 
Guenter &: Wuchterl 2003; Kenyon et al.lboosi : Joergens 2006 ). have suggested that the true 
binary fra ction of brown dwarfs may be as high as . 45 wh en selection effects are considered. 
Similarly, iPinfield et al.l (120031 ) and iGhapelle et al.l (120051 ) have inferred binary fractions of 
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0.3-0.5 for very low mass stars and brown dwarfs in two young open clusters based on the 
identification of overluminous sources in the color-magnitude diagram. 

For field brown dwarfs, particularly L/T transition objects, the identification of over- 
luminous, unresolved binaries is made difficult by uncertainties in the absolute magnitude 
relation, as w ell as intri n sic sc atter arising from age, surface gravity and radius differences. 



Nevertheless, iLiu et al.l (120061) hypothesized that a substantial fraction of L7-T5 dwarfs 
with parallax distance measurements — 10 of 15 sources (0.66lo;]^7%) — are either resolved 
or unresolved binari e s, du e to their outlying position on the Te///spectral type relation 



of iGolimowski et al.l (120041 ). Taking this suggestion at face value, and assuming that the 
current parallax sample is essentially magnitude-limited (i.e., ignoring additional selection 
effects), the fraction of known and possible binaries over this spectral type range implies 
= 0.38, with an acceptable range of 0.24-0.53 (0.21-0.59) at the 90% (95%) confidence level, 
assuminj^ P{q) oc q^. This is substantially higher than the fraction inferred from high- 
resolution imaging surveys, and suggests that up to twice as many brown dwarf binaries 
remain unresolved in these surveys. Indeed, since counting overluminous sources includes 
only those binaries with similar-mass companions, the intrinsic binary fraction may be higher 
sti ll. This result is seemingly consistent with the high brown dwarf binary fraction proposed 



bv lMaxted fc Jeffries! fl2005f ) 



Care must be taken in interpreting this value since the number of L7-T5 dwarfs with 
parallax measurements is small, uncertainties in several of the parallax measurements large, 
the underlying absolute magnitude distribution poorly constrained and hence the identifi- 
cation of overluminous sources highly uncertain. However, if the absolute magnitude scale 
is relatively fiat across the L/T transition, and the apparent fraction of binaries high, a 
measure of the fraction of overluminous sources in this spectral type regime may provide a 
more robust constraint on the true binary fraction of all brown dwarfs than high-resolution 
imaging surveys. 



6. Summary 

The primary results of this study are as follows: 

» The binary fraction excess observed amongst L/T transition objects has been suc- 
cessfully reproduced by mass function simulations incorporating the construction and 



i^For a flat mass ratio distribution, e,, = 0.42, with 0.27-0.57 (0.24-0.62) acceptable at the 90% (95%) 
confidence level. 
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classification of unresolved binary systems. It is found tliat tliis excess arises largely 
from a flattening in the M^oi/spectral type relation, which significantly reduces the 
number of single L/T transition objects in a given volume while also causing L dwarf 
+ T dwarf pairs to mimic the spectral properties of early-type T dwarfs. While the 
binary excess is particularly pronounced in magnitude-limited samples (hence its de- 
tection in current high resolution imaging studies), it is not caused by selection biases, 
but is instead intrinsic to any sample of unresolved brown dwarf systems. 

The shape of the binary fraction distribution depends weakly on the underlying ab- 
solute magnitude/spectral type relation, mass function and minimum age of the field 
population (up to 1 Gyr). The underlying binary fraction, eb produces the greatest 
effect, scaling the distribution but not significantly affecting its shape. There is a slight 
degeneracy between the influence of eb and the mass ratio distribution on the scale of 
the binary fraction distribution for a magnitude-limited sample. 

The surface density of L/T transition objects in a magnitude-limited sample is highly 
sensitive to the underlying absolute magnitude/spectral type relation. Number and 
surface densities also scale with the underlying mass function, but are largely insensitive 
to the minimum age of the field population or the mass ratio distribution of brown 
dwarf binaries. 

Empirical results constrain the underlying resolved binary fraction to lll3% (90% 
confidence interval), consistent with prior estimates. The true binary fraction, which 
includes sources unresolved in imaging studies, may be as high as 40% based on th e 



fraction of apparently overluminous L7-T5 dwarfs as proposed by iLiu et al.l (120061 ). 
However, this result requires a more robust assessment of the absolute magnitude 
scale, and understanding a selection effects in the current parallax sample, and more 
(and improved) parallax measurements of L/T transition objects. 



The simulations presented here provide new insight into the relationship between the 
fundamental properties of brown dwarfs (luminosities, mass function, intrinsic multiplicity) 
and the observed density, spectroscopic properties and multiplicity of L/T transition dwarfs. 
What they do not reveal is the physical mechanism that drives this transition. The dramatic 
change in spectral properties over a short evolutionary period is strong evidence that nonequi- 
librium dispersion of photospheric condensates must play an important, if not predominant, 
role. As new theoretical models address dynamical atmospheric effects that may be in- 
volved in condensate cloud evolution, continued empirical characterization of these sources 
— including improved multiplicity statistics through high-resolution imaging, spectroscopic 
monitoring, and parallax measurements — should be a priority. Fortunately, the apparently 
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high binary fraction of L/T transition objects implies that this remarkable evolutionary 
phase of brown dwarfs can be studied quite thoroughly, through resolved photometric and 
spectroscopic studies and long-term astrometric monitoring to measure component masses. 
These enigmatic systems may ultimately provide the most stringent empirical constraints on 
the properties of brown dwarfs in general. 
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Table 1. Empirical Binary Fraction Data 
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# Binary 
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A 1 A /I 
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''Confidence intervals based 


on the 
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distribution 



Note. 



Data compiled from hig h reso lut ion imag- 



ing surveys bv iMartin. Brandner fc Basril (119991): 



Reid et al 



( 2001bL liooil): IBouv et all J2003h: iBurpsser et ^020031. l2006al) 



Close et all (l2003l): iGizis et al. ~(l2003l) : iLiu &: LeggettI (l2005h : and 
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Table 2. Input Parameters for Monte Carlo Simulations. 



Parameter Form and Details 



Mass Function 

Age Distribution 
Evolutionary Model 



*(M) cx M-", with a = 0.0, 0.5^ , 1.0, 1.5 



^(M) oc e 



(log M-log A-lc)^ 



with Mc = O.IMq, cr = 0.627'' 



P{t) oc constant, with t e {Tq, 10 Gyr}, Tq = O.P , 0.5, 1.0 Gyr 
Baraffe et al. (2003)'' ; Burrows et al. ('20011 

0, 0.05, O.P , 0.14, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.7 



Mass Ratio Distribution P{q) oc , with F = 4'' 

P{q) oc constant 



Afboi/SpT Relation 
M^/SpT Relation 

Mj/SpT Relation 



This paper'' (see Table [3|) 

This paper** (see Table [3]) 
Liu et al. (2006) 



Liu et al. (2006 



^Parameter used in baseline simulation. 



'^Parameters from lChabrien ( 2002f ). 



Table 3. Polynomial Fits to Spectral Type Relations 



Coefficients 

Parameter Range cq ci C2 ca C4 C5 eg c^ eg cr 





LO- 
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1.374el 


1.903e-l 


1.731e-2 


7.400e-3 


-1.751e-3 
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-2.322e-6 


0.22 mag 


Mk(SpT) 
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-T8 


1.045el 


2.322e-l 
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-6.985e-6 1.051e-7 0.26 mag 


SpT(H20-J) 
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0.8 SpT 


SpT(CH4-J) 


TO- 


-T8 


2.098el 


-1.978el 


2.527el 


-3.221el 


9.087e-l 
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0.3 SpT 


SpT(CH4-K) 


LO- 


-T7 


1.885el 


-2.246el 


2.534el 


-4.734e0 


-1.259el 






1.1 SpT 



Note. — Polynomials relations are defined as f{x) = J^i^^i^''- SpT is the numerical spectral type defined as SpT(LO) = 0, SpT(L5) = 5, 
SpT(TO) = 10, etc. 
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Table 4. Spectral Templates 



Spectral Type 2MASS Photometry 

Name Opt NIR J J - Ks tt Ref 

(mag) (mag) (marcs) 

(1) (2) (3) (4) (5) (6) (7) 



2MASS J03454316+2540233 




Llitl 


14.00±0.03 


1.33±0.04 


37.1±0.5 


1,2,3 


2MASS J14392836+1929149 


Ll^ 




12.76±0.02 


1.21±0.03 


69.6±0.5 


2,4 


SSSPM J0829-1309 


L2 




12.80±0.03 


1.51±0.04 




5,6 


2MASS J15065441+1321060 


L3 




13.37±0.02 


1.62±0.03 




7 


2MASS J11040127+1959217 


L4'» 




14.38±0.03 


1.43±0.04 




8 


GJ lOOlBC'" 


L5 


L4.5 


13.11±0.02 


1.71±0.04 


77±4 


3,9,10,11 


2MASS J04390101-2353083 


L6<: 




14.41±0.03 


1.59±0.04 




8 


2MASS J09153413+0422045'' 


L7 




14.55±0.03 


1.54±0.05 




12 


DENIS J025503.3-470049 


L8 




13.25±0.03 


1.69±0.04 




13,14 


2MASS J03105986+1648155'' 


L8 


L9 


16.03±0.08 


1.71±0.11 




10,15 


SDSS J042348.57-041403.5'' 


L7.5 


rpQd 


14.47±0.03 


1.54±0.04 


65.9±1.7 


16,17 


SDSS J015141. 69+124429. 6 




Tld 


16.57±0.13 


1.38±0.23 


47±3 


15,16 


SDSS J125453.90-012247.4 


T2 




14.89±0.04 


1.05±0.06 


73.2±1.9 


16,17,18 


2MASS J12095613-1004008 




rpgd 


15.91±0.07 


0.85±0.16 




19 


2MASS J22541892+3123498 




rp4d 


15.26±0.05 


0.36±0.15 




20 


2MASS J15031961+2525196 




T5<i 


13.94±0.02 


-0.03±0.06 




21 


SDSS J162414.37+002915.6 




rpgd 


15.49±0.05 


<-0.02 


90.9±1.2 


22,23 


2MASS J07271824+1710012 


T8 




15.60±0.06 


0.04±0.20 


110±2 


17,20 


2MASS J04151954-0935066 


T8 


rpgd 


15.70±0.06 


0.27±0.21 


174±3 


17,20 



''L dwarf optical spectral standard from lKirkpatrick et al. 

^ Known binary jGoli mowski et alj l2004l : iBurgasser et al.l |2005|: iReid et alj bood: 
IStumpf. Brandner fc HennindbOOBh . 

ICruz et all ||2003| ') classify this source L6.5; however, it is adopted as an L6 spectral template here. 

'^Primary or secondary T dwarf near-infrared spectral standard from IBurgasser et al. I ll2006al') . 



Skrutski j lll997l'): (2) 'Dahn et all 1I2OO2I): (3) 



a J Scholz fc Meusirig er (2002): (6) |Lodieu et al 
(9l lGoIdman et al.l lli999i): (lOllKirkpatrick et al 



References . — (1) |Kirkpatrick, Beichman, 
iKnapp et al] ll2004h: (4) iKiA pat ricl^eXah ( 1999 
2005 ): f7') lGizis et al.l ll20|)ol): f ^Sl lCruz et al.(l200: 

2003); (ll) lHenrv et al.l l l200g) : (12) I. N. Rei d et al., in preparati on; C 13 ) iMartm et 
J. D. Kirkpatrick et al., in preparation; fl5) [ Geballe et al.l ll20 02h; 



Leggett et al. 



IVrba et al.l " ll2004h; (18) [S urgasser et al.l ll2003al'l: (19) iBurrasser et , ._ 

20021) ; (21) lBurgasser"et al. 12003b.) : (22) IStrauss et al.l (Il999l) : (23) lTinnev. Burgasser. 



112004 



(20) 



Burgasser etal 
Kirkpatrick 



2003D 
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Table 5. Number Densities as a Function of Systemic Spectral Type 





't(M) oc 










lognormal 


M-0.5 


M-0.5 




P{q) (X 


















= 


0.1 


0.1 


0.1 


0.1 


0.1 


0.3 


0.5 


LO-Ll 




0.555 


0.637 


0.712 


0.812 


0.692 


0.503 


0.396 


L1-L2 




0.755 


0.861 


1.01 


1.16 


0.973 


0.869 


0.905 


L2-L3 




0.848 


1.00 


1.16 


1.37 


1.12 


1.01 


1.01 


L3-L4 




0.693 


0.834 


1.01 


1.25 


0.962 


0.787 


0.756 


L4-L5 




0.672 


0.819 


1.03 


1.33 


0.964 


0.767 


0.707 


L5-L6 




0.688 


0.869 


1.14 


1.52 


1.03 


0.911 


0.945 


L6-L7 




0.618 


0.785 


1.04 


1.39 


0.914 


0.764 


0.746 


L7-L8 




0.502 


0.644 


0.858 


1.18 


0.772 


0.614 


0.579 


L8-L9 




0.351 


0.462 


0.607 


0.856 


0.555 


0.428 


0.393 


L9-T0 




0.237 


0.308 


0.421 


0.597 


0.366 


0.328 


0.342 


TO-Tl 




0.165 


0.220 


0.294 


0.412 


0.261 


0.251 


0.297 


T1-T2 




0.181 


0.241 


0.326 


0.442 


0.284 


0.300 


0.362 


T2-T3 




0.304 


0.406 


0.549 


0.780 


0.475 


0.452 


0.503 


T3-T4 




0.605 


0.788 


1.08 


1.54 


0.949 


0.845 


0.931 


T4-T5 




1.06 


1.42 


1.94 


2.76 


1.66 


1.40 


1.41 


T5-T6 




1.71 


2.38 


3.40 


5.06 


2.88 


2.60 


2.83 


T6-T7 




2.39 


3.45 


5.09 


7.79 


4.18 


3.23 


3.01 


T7-T8 




3.14 


4.82 


7.63 


12.4 


5.94 


4.82 


4.86 



Note. — Number densities in units of 1 0~^ pc~^ SpT~^. Results listed here are based on 
simulations using the evolutionary models of lBaraffe et al.l l|2003h , a flat age distribution spanning 
0.01-10 Gyr, the M^/spectral type relation defined in this study (§ 3.5.1), and parameters as 
specified in the header. Values are normalized to a number density of 0.00 37 pc~'^ over the mass 
range of 0.09-0.1 Mq, derived from the low-mass stellar mass function of iReid et al] l|l999l') . 



Table 6. Surface Densities for a Magnitude-Limited Sample as a Function of Systemic Spectral Type 





"^(M) oc 










lognormal 




^^-0.5 










yi-0.5 




Af(SpT):=' 


TP-K 


TP-ZC 


TP-/f 


TP-ZC 


TP-K 


L06-A" 




L06-J 


L06p-J 


TP-X 


TP-K 


TP-K 


SpT 


P{q) oc 








q* 


<?" 




q* 




q" 


1 


q* 






^6 = 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.3 


0.5 


T n T 1 




C 0/1 


O f^9 


iU.D 


19 1 


1 n Q 


fi (^9 


19/1 


1 fiT 
l.o ( 


9 /10 

z,4y 


y .00 


7 09 

/ .yz 


D.oy 


T 1 TO 




o.Do 


Q QK 

y.yo 


1 1 A 
i±.o 


1 1 


119 


1 1 

11. 


1 1 n 
11. u 


9 nn 


1 0/1 

i.y4 


o.yo 


19 9 
IZ.Z 


1 n 
10. U 


T 9 T 9 




D.Oo 


7 an 


Q QQ 


1 fl 7 


o. / U 


Q n9 


« fit; 


1 1 Q 

1. lo 


1 n9 

l.UZ 


7/11 
/ .41 


o n9 
y.uz 


1 n /I 

1U.4 


J J A 

LiO—Li^i 




O.Do 


4.00 


9fi 


(\ AR 


09 


o. ( ( 


Q /I 1 
0.41 


n /I Qn 


n /iriQ 

U.4Uo 


A 97 
4.Z / 


4. (0 


O.OO 


J A J 




9 91 


9 fiH 
Z.Oo 


O.O / 


,1 "in 

4.0U 


O. lO 


9 09 
Z.UZ 


9 (It^ 
Z.UD 


u.zuo 


n 1 Qf^ 
u. lyo 


9 RQ 
Z.Do 


Q nn 

O.UU 


Q 9Q 
o.zy 


J ^ J a 

IjO— -LiD 




1 /I V 




9 /in 




9 1 Q 

z. iy 


1 QQ 


1 /I C 


n 1 9c; 


A 1 9fi 
U.IZD 


1 i^9 
1 .OZ 


9 

Z.OO 


Q 9Q 

o.Zo 


L6-L7 




0.736 


0.932 


1.22 


1.63 


1.08 


0.818 


0.806 


0.0697 


0.0678 


0.885 


1.14 


1.34 


L7-L8 




0.398 


0.500 


0.663 


0.894 


0.600 


0.518 


0.449 


0.0497 


0.0424 


0.494 


0.639 


0.761 


L8-L9 




0.194 


0.254 


0.325 


0.464 


0.305 


0.327 


0.239 


0.0371 


0.0261 


0.245 


0.316 


0.358 


L9-T0 




0.119 


0.151 


0.205 


0.285 


0.182 


0.227 


0.130 


0.0328 


0.0201 


0.147 


0.242 


0.320 


TO-Tl 




0.0761 


0.100 


0.131 


0.186 


0.118 


0.164 


0.0777 


0.0287 


0.0163 


0.0882 


0.171 


0.242 


T1-T2 




0.0656 


0.0869 


0.117 


0.160 


0.102 


0.156 


0.0623 


0.0441 


0.0199 


0.0735 


0.150 


0.214 


T2-T3 




0.0712 


0.0952 


0.129 


0.179 


0.110 


0.171 


0.0624 


0.0738 


0.0316 


0.0857 


0.143 


0.190 


T3-T4 




0.0781 


0.102 


0.138 


0.199 


0.123 


0.179 


0.0706 


0.131 


0.0611 


0.0935 


0.134 


0.171 


T4-T5 




0.0670 


0.0897 


0.122 


0.173 


0.104 


0.142 


0.0749 


0.164 


0.0948 


0.0855 


0.102 


0.117 


T5-T6 




0.0510 


0.0711 


0.101 


0.151 


0.0859 


0.0968 


0.0705 


0.161 


0.121 


0.0644 


0.0910 


0.111 


T6-T7 




0.0307 


0.0443 


0.0652 


0.0996 


0.0536 


0.0483 


0.0463 


0.0975 


0.0926 


0.0441 


0.0462 


0.0481 


T7-T8 




0.0184 


0.0281 


0.0442 


0.0713 


0.0346 


0.0260 


0.0282 


0.0527 


0.0565 


0.0265 


0.0337 


0.0397 



00 



Note. — Surface densities in 10^^ deg"^ SpT~^. Results listed here are based on simulations using the evolutionary models of lBaraffe et al.l ||2003| ). 
a flat age distribution spanning 0.01—10 Gyr, and parameters as specified in the head er. Values are nor malized to a number density of 0.0037 pc~'^ over 
the mass range of 0.09—0.1 Mq, derived from the low-mass stellar mass function of iReid et al. 

I lll999h ; and assume a limiting magnitude of 16 in the 

appropriate photometric band. 

^Absolute magnitude/spectral type relations from (TP): this paper; (L06): iLiu et al.l ||2006| '). excluding known binaries; (L06p): iLiu et al.l ||2006|') . 
excluding known and possible binaries (see § 3.5.1). 
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Table 7. Binary Fraction as a Function of Systemic Spectral Type: Volume-Limited 





*(M) (X 


M-0.5 


M-O-5 


M-O-5 


M-O-5 


M-0.5 


M-O-5 


M-0.5 


M-0.5 


M-O-5 


M-0.5 




M(SpT):'' 


TP-K 


TP-K 


TP-K 


TP-K 


TP-K 


L06-K 


LOGp-K 


L06-J 


L06p-J 


TP-K 


SpT 


P{q) oc 




















1 




^6 = 


0.1 


0.2 


0.3 


0.4 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


LO-Ll 




0.0168 


0.0344 


0.0635 


0.0888 


0.127 


0.0122 


0.0239 


0.0178 


0.0272 


0.0455 


L1-L2 




0.111 


0.218 


0.318 


0.424 


0.526 


0.109 


0.116 


0.113 


0.117 


0.0967 


L2-L3 




0.108 


0.215 


0.311 


0.408 


0.522 


0.105 


0.101 


0.104 


0.104 


0.109 


L3-L4 




0.0788 


0.162 


0.247 


0.333 


0.435 


0.0758 


0.0726 


0.0759 


0.0780 


0.0857 


L4-L5 




0.0678 


0.137 


0.219 


0.304 


0.391 


0.0627 


0.0648 


0.0598 


0.0574 


0.0782 


L5-L6 




0.119 


0.229 


0.339 


0.446 


0.547 


0.110 


0.114 


0.106 


0.112 


0.0910 


L6-L7 




0.0911 


0.178 


0.275 


0.373 


0.467 


0.0813 


0.0919 


0.0719 


0.0763 


0.0932 


L7-L8 




0.0755 


0.154 


0.245 


0.331 


0.422 


0.0667 


0.0799 


0.0647 


0.0717 


0.100 


L8-L9 




0.0666 


0.146 


0.226 


0.309 


0.391 


0.0627 


0.0775 


0.0649 


0.0653 


0.0873 


L9-T0 




0.114 


0.228 


0.357 


0.455 


0.563 


0.115 


0.123 


0.107 


0.119 


0.128 


TO-Tl 




0.180 


0.326 


0.459 


0.573 


0.651 


0.204 


0.176 


0.184 


0.192 


0.168 


T1-T2 




0.217 


0.393 


0.521 


0.643 


0.720 


0.272 


0.208 


0.281 


0.249 


0.186 


T2-T3 




0.159 


0.294 


0.416 


0.521 


0.621 


0.187 


0.141 


0.200 


0.162 


0.138 


T3-T4 




0.139 


0.255 


0.379 


0.489 


0.587 


0.141 


0.125 


0.147 


0.133 


0.113 


T4-T5 




0.0997 


0.196 


0.297 


0.398 


0.497 


0.0983 


0.106 


0.0955 


0.101 


0.104 


T5-T6 




0.142 


0.273 


0.394 


0.502 


0.597 


0.141 


0.141 


0.141 


0.142 


0.112 


T6-T7 




0.0718 


0.147 


0.228 


0.316 


0.409 


0.0720 


0.0719 


0.0727 


0.0727 


0.0931 


T7-T8 




0.0975 


0.201 


0.298 


0.399 


0.501 


0.0995 


0.102 


0.101 


0.101 


0.101 



Note. — Results based on simulations using the evolutionary models of lBaraffe et al.l l|2003l ). a flat age distribution spanning 
0.01—10 Gyr, and parameters as indicated in the header. 

'^Absolute ma gnitude/spectral type relations from (TP): this paper; {L06): iLiu et al] ||2006| ). excluding known binaries; (L06p): 
iLiu et al] ||2006|) . excluding known and possible binaries (see § 3.5.1). 
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Table 8. Binary Fraction as a Function of Systemic Spectral Type: Magnitude-Limited 





*(M) (X 


M-0.5 


M-O-5 


M-O-5 


M-O-5 


M 


-0.5 


M-O-5 


M-0.5 


M-0.5 


M-O-5 


M-0.5 




M(SpT):'' 


TP-K 


TP-K 


TP-K 


TP-K 


TP-K 


L06-K 


LOGp-K 


L06-J 


L06p-J 


TP-K 


SpT 


P{q) oc 




















1 




^6 = 


0.1 


0.2 


0.3 


0.4 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


LO-Ll 




0.0370 


0.0737 


0.131 


0.181 





244 


0.0267 


0.0461 


0.0394 


0.0507 


0.0498 


L1-L2 




0.214 


0.378 


0.506 


0.615 





707 


0.193 


0.220 


0.198 


0.223 


0.135 


L2-L3 




0.179 


0.331 


0.449 


0.556 





665 


0.172 


0.171 


0.184 


0.188 


0.145 


L3-L4 




0.148 


0.286 


0.399 


0.506 





614 


0.158 


0.147 


0.184 


0.179 


0.122 


L4-L5 




0.150 


0.278 


0.406 


0.515 





610 


0.169 


0.150 


0.229 


0.166 


0.121 


L5-L6 




0.268 


0.453 


0.586 


0.693 





771 


0.283 


0.262 


0.311 


0.293 


0.162 


L6-L7 




0.204 


0.359 


0.496 


0.604 





692 


0.215 


0.204 


0.250 


0.220 


0.169 


L7-L8 




0.214 


0.381 


0.521 


0.622 





710 


0.180 


0.197 


0.238 


0.249 


0.212 


L8-L9 




0.200 


0.372 


0.506 


0.609 





685 


0.167 


0.207 


0.225 


0.224 


0.183 


L9-T0 




0.317 


0.530 


0.671 


0.755 





824 


0.264 


0.327 


0.314 


0.362 


0.300 


TO-Tl 




0.428 


0.621 


0.746 


0.820 





863 


0.387 


0.472 


0.354 


0.452 


0.369 


T1-T2 




0.435 


0.640 


0.751 


0.831 





877 


0.453 


0.499 


0.421 


0.466 


0.356 


T2-T3 




0.330 


0.525 


0.654 


0.741 





812 


0.351 


0.365 


0.281 


0.294 


0.255 


T3-T4 




0.249 


0.413 


0.558 


0.662 





747 


0.242 


0.256 


0.248 


0.257 


0.173 


T4-T5 




0.174 


0.316 


0.441 


0.556 





649 


0.158 


0.191 


0.161 


0.187 


0.142 


T5-T6 




0.228 


0.397 


0.535 


0.640 





723 


0.206 


0.231 


0.220 


0.243 


0.149 


T6-T7 




0.119 


0.233 


0.343 


0.448 





550 


0.112 


0.119 


0.116 


0.122 


0.119 


T7-T8 




0.185 


0.346 


0.472 


0.582 





678 


0.189 


0.194 


0.189 


0.189 


0.140 



Note. — Simulations based on the evolutionary models of lBaraffe et al ] l|2003h . a flat age distribution spanning 0.01-10 Gyr, 
and parameters as indicated in the header. A limiting magnitude of 16 in the respective filter band is assumed. 

'^Absolute ma gnitude/spectral type relations from (TP): this paper; (L06): iLiu et al] ||2006| ). excluding known binaries; (L06p): 
iLiu et al] ||2006|) . excluding known and possible binaries (see § 3.5.1). 
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LO L2 L4 L6 L8 TO T2 T4 T6 T8 

Spectral Type 



Fig. 1. — Bolometic magnitude versus spectral type for 26 unresolved field sources (open 
circles) and the six components of the binaries Kelu lAB, Epsi l qn Indi Bab and SPSS 0423 - 
0414AB (filled circles ). D ata are fromlGolimowski et all (120041 ): McCaughrean et all (120041 ): 
Liu fc LeggettI (120051 ): and lBurgasser et all (|2006 J T" A sixth-order polynomial fit to the data, 



excluding the apparently overluminous sources SDSS 1254-0122 and 2MASS 0559-1404, is 
indicated by the solid line (see Table [3]). The residual scatter in the fit is 0.22 mag. 
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Fig. 2. — SpeX spectral templates for L dwarfs (left) and T dwarfs (right). Spectra are 
normalized in the 1.1-1.3 fim region and offset by constants for clarity (dotted lines). Names 
and spectral types are listed (see Table HI). 
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LO L2 L4 L6 L8 TO T2 T4 T6 T8 

Spectral Type 



Fig. 3. — Absolute MKO i^-band magnitude versus spectral type for 26 unresolved field 
sources (open circles) and the eight components of the binaries Kelu lAB, Epsilon Indi 



Geballe et al. 




^002): 


Leeeett et al. (2002): 


KnaDD et al. 


(2004): McCauehrean e 


t ah (2004): 


Liu & LeeeettI 


(2005 


); and 


Bureasser et al 


. ( 


2006a). Parallax data are from 


Dahn et al. 


(2002h: Tinnev. Bureasser, & Kirkpatrick ( 


2003): and ^ 


Vrba et al. (2004). A eiehth-order 



polynomial fit to these data is indicated by the black line (see Table [3]). M/^/spectral type 



relations from iLiu et al.l (120061 ). excluding known binaries (red dashed line) and known and 



possible binaries (blue dot-dashed line) are also shown. 



-42 - 




L4 L6 L8 TO T2 T4 T6 

Primary Spectral Type 



Fig. 4. — Hybrid binary spectra for various combinations of primary (ordinate) and secondary 
(abscissa) spectral type combinations. All spectra are normalized in the 1.1-1.3 /xm band. 
Composite spectral types, based on spectral indices as described in § 3.5.2, are indicated in 
the upper right corner of each box. 
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Fig. 5. — Spectral indices as a function of spectral type for a large sample of SpeX prism 
data. Measurements made on the spectral templates shown in Figure [2] are indicated by 
filled circles, all others by open circles. Spectral types are based on optical data for L0-L8 
dwarfs and near-infrared data for L9-T8 dwarfs. Fourth-order polynomial fits for each index 
over a prescribed range are indicated by solid lines (see Table |3]) . 
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Fig. 6. — Examination of the robustness of near- infrared classifications for L and T dwarfs 
using the spectral index relations of Table [31 Plotted in small open circles are the derived 
spectral types for 88 L0-T8 dwarfs with SpeX prism data using the spectral index relations 
of Table [3|, compared to their published optical (L0-L8) and near-infrared (L9-T8) spectral 
types. Average deviations and la scatter from perfect agreement (dotted line) in groupings of 
two spectral subtypes are indicated by filled circles with error bars. The typical disagreement 
for the entire sample is 0.6 SpT, but is greater for L0-L8 dwarfs (0.9 SpT) than L9-T8 dwarfs 
(0.3 SpT). 




Fig. 7. — Number ($(SpT); left panel) and surface density (S(SpT); right panel) distri- 
butions for single sources based on simulations using power-law mass functions with a = 
(red dashed line), 0.5 (solid black line), 1.0 (blue dot-dashed line ) and 1.5 (green triple- 
dot-dashed line), and the lognormal distribution of IChabrierl (j2002l . magenta dotted line). 
Number densities are shown in units of pc^^ SpT^^ for a volume-limited sample, surface 
densities are shown in units of deg'^ SpT^^ for a magnitude-limited sample with K < 16. 
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Fig. 8. — Number (left panel) and surface density (right panel) distributions for primaries 
(solid black lines), secondaries (blue dot-dashed hues) and systems (red dashed lines) for the 
baseline simulations. Surface densities are based on a magnitude-limited sample {K < 16). 
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Fig. 9. — Surface density distributions as a function of systemic spectral type for different 
absolute magnitude/spectral type relations. (Left panel) M/^/spectral type relations from 
this study (black line; Table [3l) and lLiu et al.l (120061 ). excluding known (red dashed line) and 
known and possible binaries (blue dot-dashed line). (Ri ght panel) Comp arison of Mj (red 
lines) and (black lines) spectral type relations from iLiu et al.l (120061 ) excluding known 
(solid line) and known and possible binaries (dashed lines). All distributions are based on a 
magnitude-limited sample with a limiting magnitude of 16 in their respective bands. 
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Fig. 10. — Binary fraction distributions for volume-limited (left) and magnitude-limited 
samples (right). In both panels, the binary fraction as a function of systemic spectral type 
is indicated in black, while the input binary fraction distribution is indicated in red {eb = 
0.1). 
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Fig. 11. — Binary fraction distributions for a magnitude-limited sample with = 0.05 (red 
dashed line), 0.1 (black line), 0.15 (blue dot-dashed line), 0.2 (green triple-dot dashed line) 
and 0.3 (magenta dotted line), resampled into spectral type bins of L0-L4, L4-L8, L8-T4, 
and T4-T8. Empirical measurements are indicated by solid circles with error bars, where 
the uncertainties corr espond to the 90% confidence intervals for a binomial distribution 
( jBurgasser et al.ll2003l ). 
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Fig. 12. — Binary frac t ion d istributions for magnitude-limited samples assuming different 
forms of the Liu et al. Mk (black lines) and Mj (red lines) spectral type relations, 

excluding known (solid lines) or known and possible binaries (dashed lines). The left panel 
shows the distributions sampled by integer spectral type; the right panel shows the same 
distributions sampled by the spectral bins listed in Table [H Empirical measurements are 
indicated in the right panel as in Figure [Til 
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Fig. 13. — Binary fraction distributions for magnitude-limited samples assuming different 
forms of the underlying mass ratio distribution. Black lines plot distributions assuming e?, 
= 0.1 and P{q) oc g^, red dashed lines plot distributions assuming eb = 0.1 and P{q) oc 
constant, and blue dot-dashed lines plot distributions assuming = 0.14 and P{q) oc 1. 
The left panel shows the binary fraction distributions sampled by integer spectral type; the 
right panel shows the same distributions sampled by the spectral bins listed in Table [H 
Empirical measurements are indicated in the right panel as in Figure [HI 
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Fig. 14. — Composition of binaries as a function of composite spectral type. (Top panel) 
Mean primary spectral type, where the vertical error bars indicate the minimum and max- 
imum primary spectral type of the binaries for a given composite spectral type. (Middle 
panel) Average difference between composite and primary spectral type. (Bottom panel) 
Percentage of binaries with primaries classified a full spectral class or earlier than the com- 
posite type. 



